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Abstract 

Test particle transport determined by the Lorentz force in turbulent mag- 
netized plasmas is studied. The time dependent diffusion coefficient, valid 
for the whole range of parameters, is obtained by developing the decorrela- 
tion trajectory method. The effects of Larmor gyration in the presence of 
trapping are analyzed and several transport regimes are evidenced. 
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1 Introduction 

Impurity control in magnetically confined plasmas is a very important issue for the 
development of fusion reactors. Impurity behaviour is a complex problem related to 
confinement and transport of the bulk ions and electrons in plasma and to plasma- 
wall interaction. A very strong experimental effort (see e.g. PP-[0]) lead to the 
conclusion that this process is far from being understood on the basis of the ex- 
isting theoretical models. In particular, the experimental results for the diffusion 
coefficients are much larger than the neoclassical prediction, especially in the edge 
plasma, showing the presence of anomalous transport. 

We analyze here the particular topic of impurity transport in turbulent plasmas 
using test particle approach. Particle motion in a stochastic potential was exten- 
sively studied in the guiding center approximation ^UJ- . It is well known since 
many years jTHj that, for slowly varying or large amplitude turbulence, the E x B 
drift determines a process of dynamical trapping of the trajectories. It consists of 
trajectory winding around the extrema of the stochastic potential and strongly in- 
fluences the transport. Important progresses in the study of this nonlinear process 
were recently obtained. New statistical methods were developed [T7j, [TH] that per- 
mitted to determine the asymptotic diffusion coefficient and also the correlation of 
the Lagrangian velocity and the time dependent (running) diffusion coefficient. It 
was shown that the trapping process determines the decrease of the diffusion coef- 
ficient and the change of its scaling in the parameters of the stochastic field. These 
methods were extended to more complex models, which consider, besides the E x B 
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drift, particle collisions |19j . an average velocity [20] or the parallel motion (21], and 
also to the study of the collisional particle diffusion in stochastic magnetic fields 
[22J. The conclusion of these studies is that the trapping combined with a decor- 
relation mechanism determines anomalous transport regimes. In these regimes the 
dependence of the diffusion coefficient D on the parameters describing the decor- 
relation is inverted due to trapping in the sense that a decrease of D appearing 
in the quasilinear regime is transformed into an "anomalous" increase of D in the 
nonlinear regime [2*3] . 

All these studies are based on the guiding center approximation for particle 
motion that considers the Larmor radius negligible. This approximation is not 
adequate for the impurity ions which can have Larmor radii comparable or larger 
than the correlation length of the turbulence and cyclotron periods comparable with 
the turnover time of the E x B motion. In these conditions the trajectories have to be 
determined from the Lorentz force. The aim of this paper is to determine the effects 
of finite Larmor radius on particle transport in a turbulent magnetized plasmas. 
In particular, we analyze the influence on the trapping process and compare the 
characteristics of the transport induced by the Lorentz force (Lorentz transport) 
with those obtained in the guiding center approximation (drift transport). The time 
dependent diffusion coefficient is obtained as a function of the turbulence and ion's 
parameters by extending the decorrelation trajectory method developed for the drift 
transport [IZ|, [23- The transport regimes for a large range of parameters are 
determined. 

The paper is organized as follows. Section 2 contains the basic equations, the 
statistical approach and the derivation of the general expression for the time de- 
pendent diffusion coefficient. The results are presented in Sections 3-6. First, in 
Section 3, a static potential is considered. We show that the trapping exists even at 
large values of the Larmor radius and that it determines a subdiffusive transport, as 
in the drift approximation, but with a time dependence of the diffusion coefficient 
strongly influenced by the gyration motion. Then we examine the dependence of the 
asymptotic diffusion coefficient on the three dimensionless parameters (see Section 2 
for their definitions) that characterize this process and show that several transport 
regimes appear. The dependence on the first parameter (the Kubo number) that 
describes the effect of the time dependence of the stochastic potential is analyzed 
in Section 4 where two regimes are evidenced. The dependence of D on the initial 
kinetic energy of the ions (the second parameter) is presented in Section 5 and in 
Section 6 the effect of the normalized cyclotron frequency that essentially describes 
the specific charge of the ions is examined. The conclusions are summarized in 
Section 7. 

2 Basic equations and statistical approach 

We consider a constant confining magnetic field directed along z axis, B = Be z (slab 
geometry) and an electrostatic turbulence represented by an electrostatic potential 
0'(x, t), where x = (x±, x%) are the Cartesian coordinates in the plane perpendicular 
to B. The motion of an ion with charge q and mass m is determined by the Lorentz 
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force: 

m ^T = Q {" W'(x, t) + u x B} (1) 

where x(i) is the ion trajectory, u(i) = dx(t)/dt is its velocity and V is the gradient 
in the (xi,X2) plane. The initial conditions are 

x(0) = 0, u(0) = u . (2) 

This equation is transformed into a system of first order equations for the position 
and the velocity of the ion 

dui q d0'(x,t) 

dx\ . . 

IF = "< (4) 

where f2 = qB/m is the cyclotron frequency and £j„ is the antisymmetric tensor 
(e 12 = — e 2 i = 1, £n = £22 = 0). Introducing the instantaneous Larmor radius 
defined by 

Uj(t) 

n 

the guiding center position £(£) = x(£) — p(t) and 0(x, t) = (f>'(x,t)/B, the system 
becomes 



PM = (5) 



dp* _ 

dt ij 



d^ + p.t) 

- + n Pj 



(7) 



The electrostatic potential 0(x, £) is a stochastic field and thus Eqs. (JBHZJ) are 
Langevin equations. The solution consists, in principle, in determining the statis- 
tical properties of the ensembles of trajectories, each one obtained by integrating 
Eqs. (jMZj) for a realization of the stochastic potential. We will determine here the 
mean square displacement and the time dependent diffusion coefficient for the guid- 
ing center trajectories £(£). These statistical quantities can also be determined for 
particle trajectories x(t) and for the Larmor radius p(t) but they are not physically 
relevant. 

The potential is considered to be a stationary and homogeneous Gaussian stochas- 
tic field, with zero average. Such a stochastic field is completely determined by the 
two-point Eulerian correlation function, E(x.,t), defined by 

£(x, t) = (0(x\ t>) 0(x- + x, t> + t)) . (8) 

The average (...) is the statistical average over the realizations of 0(x, t). The sta- 
tistical properties of the drift velocity components 



. (9) 
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are completely determined by those of the potential; they are stationary and ho- 
mogeneous Gaussian stochastic fields like 0(x, t). The two-point Eulerian correla- 
tions of the drift velocity components, (x, t) = (vf r {x%t') t^ r (x' + x, , and 
the potential-velocity correlations, E^x, t) = (</>(x', £') t>f r (x' + x, t' + £)) , are ob- 
tained using Eq. © as: 

d 2 E(x,t) 

Eij{X,t) — —e in ej m — — , (10) 

9E(x,t) 

E$i{x,t) = —Ei ( j > (x,t) = —e in — — . 

The Eulerian correlation of the drift velocity evidences three parameters: the 
amplitude V = a/ -E'ii(0, 0), the correlation time r c , which is the decay time of the 
Eulerian correlation, and the correlation length A c , which is the characteristic decay 
distance. These parameters combine in a dimensionless Kubo number 

A c Tfi 

which is the ratio of r c to the average time of flight of the particles (r ^ = A c /V) 
over the correlation length. Using these parameters of the stochastic field, Eqs.(07J) 
are written is dimensionless form as 

d^__ c <90(£ + p,t) 

dt ~ £lJ [ ] 



dt ij 



d<j>(d + p.t) n 



(13) 



where the normalization parameters are r f\ for time, A c for distances, V for the drift 
velocity and 

tt = ttT fl . (14) 

The same notations are kept for the normalized quantities. 

We note that the equation for the guiding center trajectory (fT^j) is similar with 
that obtained in the guiding center approximation, with the difference that the 
argument of the potential is the particle trajectory x(t) = £(t)+p(t) instead of £(t). 
The equation for the Larmor radius (jlHjl describes a cyclotron motion that has the 
radius and the frequency dependent of the stochastic potential. In the 2-dimensional 
case studied here, both equations are of Hamiltonian type (the two components of 
£(t) are conjugated variables as well as the two components of the Larmor radius 
Pit) 

dH(£,p) d Pl _ dH((j, p) 

dt - ij d£j ' dt 13 d pj • 1 } 

They have the same Hamiltonian function 

H(tP) = m+p) + ^(pi+pDi as) 

which is the energy of the particle. The Hamiltonian depends on £ and p and thus 
the two Hamiltonian systems (fT2|) and (fT3*|) are coupled. For each system the other 
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variable introduces a time dependence in H (£, p) which perturbs the regular motion 
that is obtained in the absence of interaction. The perturbation can be very strong 
leading to a chaotic motion of the guiding centers. 

Particle motion is thus determined by three dimensionless parameters: K, ~p and 
Q. The first one, the Kubo number K, does not appear in the equations, but only in 
the statistical description of the stochastic potential. It describes the effect of time 
variation of the stochastic potential. The second parameter p is the initial Larmor 
radius normalized with the correlation length and appears in the initial condition 
(p|) . which is written as 

Pl (0) = pcos(a), p 2 (0) = psin(a), £(0) = -p(0) (17) 

where p = = |u 1 /VVt and a determines the orientation of the initial velocity 
(the angle between u and the x\ axis is 7r/2 — a), p is related to the initial kinetic 
energy of the particles. The third parameter Q defined in Eq. (jl4)l is the cyclotron 
frequency normalized with and describes the relative importance of the cyclotron 
and drift motion (second and respectively first term in I13j)) in the evolution of the 
Larmor radius. 

Starting from the statistical description of the stochastic potential, we will de- 
termine the correlation of the Lagrangian drift velocity, defined by: 

= <«f [x(0),0])v*[x(f),f]). (18) 

The mean square displacement of the guiding center and its time dependent diffusion 
coefficient are integrals of this function: 

<C 2 (t)>=2 f dTLn(r)(t-T), (19) 
Jo 

Di(t)= [ drL«(r), (20) 
Jo 

provided that the process is stationary 12 ")j . 

The guiding center approximation obtained by taking p = in Eq. (fT2"|) was re- 
cently studied by developing a semi-analytical approach, the decorrelation trajectory 
method, [23]. Using this approach an important progress was obtained in the 
understanding of the intrinsic trapping process specific to the E x B drift. We 
present here a generalization of the decorrelation trajectory method that applies to 
the Lorentz transport described by Eqs. (| 121113)) . 

The Langevin equations (jl2H13p for given values of the parameters Q, ~p and 
K is studied in subensembles (S) of realizations of the stochastic field, which are 
determined by given values of the potential and of the drift velocity in the starting 
point of the trajectories: 

0(0,0)=/, v^(0,0) = v°. (21) 

All the trajectories contained in a subensemble have the same initial energy. The 
stochastic (Eulerian) potential and drift velocity in a subensemble (S) defined by 
condition ()21j) are Gaussian fields but non-stationary and non-homogeneous, with 
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space and time dependent averages. These averages depend on the parameters of 
the subensemble and are determined by the Eulerian correlation of the potential 



$(x,i;S) = (0(x, t)) ? = u ^Hf + v" LV ) ' + ^ 2 U ' , (22) 
v , , ; ^ » v/s ^£(0,0) ^11(0,0) 2 £ 22 (0,0)' v ; 

V^(x,t; 5) = < V * [x, t]) s = ~e l3 m ^ S) (23) 

where (■■■)$ is the statistical average aver the realizations that belong to (S). They 
are equal to the corresponding imposed condition (J2~T| in x = and t — and decay 
to zero at large distance and/or time. The existence of an average Eulerian drift ve- 
locity in the subensemble determines an average motion, i.e. an average Lagrangian 
drift velocity V L (t; S) = (w dr [x(i),£]^ . The correlation of the Lagrangian drift 
velocity for the whole ensemble of realizations (fTSj) can be written as 

Lij(t) = J Jd<P° dv° P 1 ( ( p , v°; 0, 0) v^ L (t; S) (24) 

where -Pi(0°, v°; 0, 0) is the probability that a realization belongs to the subensemble 
(S). The average Lagrangian drift velocity V L (t; S) is determined using an approx- 
imation that essentially consists in neglecting the fluctuations of the trajectories in 
(S). This approximation is validated in ^Hj where the fluctuations of the trajectories 
in (S) are taken into account in a more complicated and precise method. It is shown 
that they lead to a weak modification of the diffusion coefficients D(t), although 
they strongly change V L (t; S). Introducing the average guiding center trajectory in 
(S), E(£; S) = (£(t)) s , and the average Larmor radius in (S), H(t; S) = (p(t)) s , the 
equations of motion can be averaged over the realizations in (S) in this approxima- 
tion and yield 

d~ t d^(3 + U,t;S) 

dt ~ £i3 9E< ' 



djk _ 

dt 



(26) 



The initial conditions for the two components of the subensemble average trajectory 
are obtained from Eq. (fTTj) 

rii(0) = pcos(a), n 2 (0) = psin(a), 3(0) = -11(0). (27) 

Since the orientation of the initial velocity Uo is arbitrary, we will consider in each 
realization in (S) initial conditions with uniform distribution of a over the interval 
[0,2vr). 

This approximation ensures the conservation of the subensemble average energy 
of the particles. 

Considering for simplicity an isotropic stochastic potential with the Eulerian cor- 
relation depending |x| , a diagonal correlation tensor is obtained for the Lagrangian 
drift velocity, = SijL(t), and the following expressions for the time dependent 

(running) diffusive coefficient: 

D(t) = D B F(t), (28) 
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Hi) = v- 



2 dF(t) 
dt 



(29) 



where 



Fit) 



2(2tt) 



1 



3/2 




da Si(t; S) 



(30) 



and 5i(t; 5") is the component of the solution of Eqs. (j25H2l)|) along the initial average 
drift velocity v° and v — |v°| . 

We have thus determined the correlation of the Lagrangian drift velocity (for 
ions with mass m, charge q and given initial kinetic energy) corresponding to given 
Eulerian correlation £"(x, t) of the stochastic potential. Explicit results for L[t) and 
D[t) are obtained by effectively calculating the average trajectories in (S), solutions 
of Eqs. (|25H27|) . and the weighted average (JHOJ). This procedure appears to be 
very similar with a direct numerical study of the simulated trajectories. There are 
however essential differences. The average trajectories are obtained for a rather 
smooth and simple potential and the number of trajectories is much smaller than 
in the numerical study due to the weighting factor determined analytically. This 
reduced very much the calculation time, such that it can be performed on PC. A 
compuder code is developed for explicit calculation of D{t) for given values of the 
parameters K, ~p, Q and prescribed Eulerian correlation of the potential. 

The results presented in next sections are for -E^x, t) = 1/(1 + x 2 /2) exp(— t/K). 
As shown in [23] the shape of the Eulerian correlation of the potential determines the 
strength of the trapping represented by the exponent of the time decay of the diffu- 
sion coefficient in the static case. However the general behavior of the decorrelation 
trajectories and of L(t) and D(t) are the same for all correlations. 

3 Transport by Lorentz force 

We consider here a static potential (K, t c — > oo) and compare the results with those 
obtained in the guiding center approximation (drift transport) jTJj, [231 ■ The aim 
is to identify the effects of the Larmor radius. 

In the frame of the decorrelation trajectory method, the difference between the 
Lorentz and the drift transport consists in the equations for the decorrelation tra- 
jectories. The trajectories obtained from Eqs. (j25H26|) are much complicated than in 
the guiding center approximation. The important simplification introduced in the 
latter case that actually reduces the number of parameters at one and eliminates 
the time dependence of the stochastic potential (see [T^j, [231) cannot be applied in 
the present case. The trajectories effectively depend on the six parameters a, u, 0°, 
Q, p and K. In the static case there are closed periodic trajectories of the guiding 
center £(£), even at large Larmor radii which shows that trapping exists. The initial 
drift velocity does not influence only the period but also the size and the shape of 
the paths. The orientation of the initial velocity of the particle, a, strongly influ- 
ences the trajectories obtained with the initial condition f!27|) because the average 
produced by the gyration is different for different values of a. 
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Figure 1: The contour plot of the function /(</> ,t) for the drift (a) and Lorentz (b) 

transport. 

The function F(t) in Eq. 1)30)1 is obtained by summing the contribution, f(<fi°, t), 
of all the decorrelation trajectories that start from a point where the potential is 
<fi°. This function t) gives details of the diffusion process showing the trapping. 
The contour plot of this function is represented in Fig. l.a. for the drift transport 
and in Fig. l.b. for the Lorentz transport. One can see that in the first case the 
contributions of the large values of are eliminated progressively due to trapping 
and f{4>° 1 t) shrinks continuously as time increase such that thinner and thinner 
intervals of 0° centered around zero contribute to D{t). The Lorentz transport is 
characterized by a completely different pattern for this function (Fig. l.b). At 
small time it has a Gaussian shape that increase up to a maximum. This increase 
is not smooth but performed in steps. At later times, / becomes unsymmetrical 
and has a minimum at (jP = 0. It decays continuously in time maintaining a large 
range of (fP . Thus, the trapping process is completely different for the Lorentz 
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transport. At each value of in a large interval, the cyclotron motion determines 
large, open trajectories as well as small, closed ones, depending on the values of a 
and v. The contributions of small trajectories are progressively eliminated by mixing 
determining the decay of the function f((j)°,t). 

The time dependent diffusion coefficient is presented in Fig. 2 for the Lorentz 
transport (continuous line) and for the drift transport (dashed line). The time 
dependence of the Lorentz diffusion coefficient is rather complex and a strong influ- 
ence of the Larmor radius can be observed. At small time the diffusion coefficient 
increases nonuniformly, in steps. Averaging these steps a linear time dependence 
can be observed, similar with that obtained in the drift transport. This behavior 
extends to times much longer than the flight time. The maximum of D(t) is at 
about 7tj[. At later times a decay of D(t) appears with a time dependence that is 
approximately the same as in the drift case, but with D(t) larger with a factor of 
about 2. Thus the transport in the static case is subdiffusive. 




-2 -1 1 2 

10 10 10 10 . 10 



Figure 2: The time dependent diffusion coefficient for the Lorentz transport (blue 
line) and for the drift transport (dashed line) for Q = 10, p = 1, K = oo. 

The correlation of the Lagrangian drift velocity is presented in Fig. 3 for the 
Lorentz transport (continuous line) compared with drift approximation (dashed 
line). It decays very fast (in a time much smaller than Tfi) and then it presents 
a series of peaks with decreasing amplitude and eventually has a negative tail. The 
peaks appear around multiples of the cyclotron gyration period T = 2n/fl. 
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Figure 3: The correlation of the Lagrangian drift velocity for the Lorentz transport 
(blue line) and for the drift transport (dashed line) for Q — 10, p = 1, K — oo. 

A clear story of the physical process can be deduced from the time evolution 
of D(t) and L(t). Starting from t — 0, at very small time (t <C T) D(t) is equal 
with the drift diffusion coefficient. Then, the cyclotron motion with a large Larmor 
radius ( p — 1 in Figs. 2, 3) averages the stochastic field along the trajectory. 
Consequently the guiding center has a very small displacement and D(t) is much 
reduced compared to the drift case. After a period the trajectories come back near 
the initial position (all in phase because T is a constant), a coherent motion of 
the guiding centers appears during the passage of the particles and a step in D{t) 
is produced. Thus the evolution of the guiding centers is determined mainly by 
short coherent kicks appearing with period T. Their displacement is thus slower. 
Consequently, the trapping appears at later time and the linear increase of D(t) 
extends to longer times leading to values of D(t) that are higher than for the drift 
diffusion. When the displacements of the guiding centers increase the coherence of 
the periodic kicks is progressively lost and the peaks of L{t) becomes smaller and 
thicker. 

The diffusion coefficient obtained for the Lorentz transport is a function of the 
dimensionless parameters K, ~p and Q. They contain the physical parameters of the 
stochastic field (amplitude, correlation length, correlation time) and of the impurity 
ions (mass, charge, kinetic energy). The dependence of the diffusion coefficient on 
these parameters is analyzed in the next Sections. 

4 K dependence 

We consider here time dependent stochastic potentials and determine the depen- 
dence of the diffusion coefficient on the parameter K defined in Eq. flllj) . 

In the case of drift transport, a change of variable can be done in the equation for 
the decorrelation trajectories in order to introduce the time factor of the Eulerian 
correlation of the potential in the time variable. The diffusion coefficient for the time 
dependent potential is so determined from D(t) obtained for the static potential with 
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the same space dependence in the Eulerian correlation [2*lj . The equations for the 
decorrelation trajectories (j25H26|) obtained for the Lorentz transport do not have 
this property: due to the cyclotron motion (second term in Eq. (|26|)). the time 
factor in the average potential $ cannot be introduced in the time variable. Thus 
the K dependence of the diffusion coefficient must be determined by performing the 
calculations of D(t) for each value of K. 





10" 2 10"' 10° 10' t/t, 10 2 



b 

Figure 4: The evolution of the diffusion coefficient for the Lorentz transport in a 
time dependent stochastic potentia (red line) compared to the results obtained in a 
static potential for Lorentz transport (blue line) and for the drift transport 

(dashed line). 

However, we have shown in ^H] , using general considerations based on the shape 
of the correlation of the Lagrangian velocity, that at large K the asymptotic diffusion 
coefficient can be approximated by the diffusion coefficient determined for the static 
potential at t — r c , (or in dimensionless units at t — K). These considerations are not 
dependent on the specific type of motion or on the statistical method used to obtain 
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the Lagrangian correlation. We examine here the accuracy of this approximation 
for the Lorentz transport. 

Typical examples for the evolution of the diffusion coefficient in time dependent 
stochastic fields are presented in Fig. 4 (continuous line) compared to D{t) obtained 
in the static potential for Lorentz (dotted line) and drift (dashed line) transport. 
One can see that the diffusion coefficient saturates showing that the transport is 
diffusive in time dependent stochastic potentials. A large K case is considered 
in Fig. 4. a which shows that the above approximation is rather accurate. For 
smaller K values that are not situated on the tail of the Lagrangian correlation, the 
demonstration presented in ^H] does not apply. However, as seen in Fig. 4. b for 
K = 2 the above approximation is valid and so was in many other cases we have 
considered. Thus the asymptotic diffusion coefficient in a time dependent stochastic 
potential with Kubo number K is 

D(oo\K,p,0) = D(K\oo,p,0) (31) 

where D(t\K, ~p, Q) is the time dependent diffusion coefficient obtained for the param- 
eters K, p, Q, D(oo\K,p, Q) is its asymptotic value and D(t\oo,p, Q) is the diffusion 
coefficient obtained in the static potential. 

5 p dependence 

The parameter p = |p(0)| /A c essentially describes the effect of the initial kinetic 
energy of the ions on their Lorentz diffusion. It does not appear in the drift transport, 
which is determined only by the Kubo number K. 

The dependence of the asymptotic diffusion coefficient on the Kubo number for 
several values of p is presented in Fig. 5. One can see that the Larmor radius 
produces observable effects even for rather small values (at p = 0.1). As expected, 
the effect strongly increase with the increase of p. This modification of the diffusion 
coefficient due to Larmor radius is complex and it may consist of a strong decrease 
as well as of a strong increase, depending on the conditions. Thus, the general idea 
that the effective diffusion is reduced due to the cyclotron motion which averages 
the stochastic potential, is not always true. 
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Figure 5: The asymptotic diffusion coefficient for the Lorentz transport as a 
function of K for several values of p that label the curves and for Q = 10. The 
result obtained in the guiding center approximation is also represented (dashed 

line) . 

At small Kubo numbers the diffusion coefficient is much smaller than in the drift 
approximation. It increases in steps that appear, for all values of p, at values of K 
which are multiples of the cyclotron period T = 2n/Q. Apart from these steps that 
are attenuated at larger K, there is a global increase with K as D{po\K, p, fl) ~ 
DbK = (\1/t c )K 2 in this regime. This is similar with the quasilinear regime of 
the drift transport and corresponds to initial ballistic motion of the guiding centers. 
This regime extends to K > 1, up values that increase with the increase of p. The 
diffusion coefficient at this value of K has a value much larger than for the drift 
transport. 

At larger values of K, the trapping becomes effective and the diffusion coefficient 
has approximately the same K dependence as in the drift transport. The effect of 
the Larmor radius consists in an amplification factor in the diffusion coefficient that 
is independent of K in this regime. It increases with the increase of ~p. 

6 Q dependence 

The parameter Q = Qrfi describes the effect of the specific charge of the ions q/m 
on their Lorentz diffusion. It determines the moments of the steps appearing in the 
K dependence of the diffusion coefficient. Apart this, there is no strong influence 
of this parameter when Q ^> 1 (see Fig. 6). For small Q the trajectories become 
chaotic and the diffusion coefficient has an irregular dependence on the parameters. 
Note that the chaotic variations seen in this Figure at large K for Q=l are not 
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calculation errors (they are not changed by increasing the number of calculated 
trajectories, thus the accuracy). 




10" 2 10" 1 10° 10 1 



Figure 6: The asymptotic diffusion coefficient for the Lorentz transport as a 
function of K for the values of Q that label the curves and for p — 1; the result 
obtained in the guiding center approximation is also represented (dashed line). 

7 Conclusions 

We have studied the impurity ion transport produced by the Lorentz force in a tur- 
bulent magnetized plasma. Expressions for the time dependent diffusion coefficient 
and for the correlation of the Lagrangian drift velocities are obtained in terms of 
a class of smooth, deterministic trajectories by developing a generalization of the 
decorrelation trajectory method. This statistical approach is compatible with the 
invariance of particle energy. 

We have shown that the Larmor radius has a strong effect on impurity ion trans- 
port in turbulent plasmas. The generally accepted idea that the effective diffusion is 
reduced due to the cyclotron motion which averages the stochastic potential, is not 
always true. The cyclotron motion can also determine the build up of correlation of 
the Lagrangian drift velocity by bringing the particles back in the correlated zone of 
the stochastic potential. The correlation L(t) shows a series of periodic peaks, which 
lead to incresed diffusion coefficients in slowly varying potentials. Consequently, at 
given Larmor radius, the transport can be reduced or increased, depending essen- 
tially on the value of the Kubo number. 
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